Boundary layer parameterization scheme is one of the most important sources of numerical model predict error. Selecting the suitable planetary boundary layer parameter scheme for the simulation of the Bohai costal area typhoon rainstorm can provide reference for the subsequent business application and scientific research in this area. Eight planetary boundary layer (PBL) parameterization schemes (ACM2, BouLac, GBM, MYJ, MYNN, QNSE, UW, YSU) in the Weather Research and Forecast Model (WRFV4.3) are used to simulate the torrential rain caused by northward Typhoon In-Fa (2106). The discrepancies of the simulated torrential rain based on the different schemes are compared, and the simulated PBL thermo-dynamical structure are examined based on ERA5 data. The results are as follows: Firstly, the simulated typhoon rainfall (including 24 h accumulated rainfall, the maximum and position of the accumulated rainfall, equitable threat scores, the hourly maximum rain and the rainfall areas of 10.0 mm and 20.0 mm precipitation per hour) during the typhoon northward period present significant difference, and the simulated track difference mainly occurred in the middle-late simulation period. Secondly, the local closure scheme-BouLac test has the best equitable threat scores for the 24 h accumulated rainfall above 10.0 mm, while the non-local closure scheme-ACM2 simulated the accumulated rainfall above 25.0, 50.0 and 100.0 mm best, which has the highest equitable threat score compared to the other PBL parameterization schemes. In addition, simulation results of the ACM2 scheme show the maximum accumulated rainfall, 24 h accumulated rainfall averaged on a region and mean hourly maximum rainfall are closest to ERA5 data. Therefore, the ACM2 scheme is best suitable for the torrential rain simulation of Typhoon In-Fa (2106) under the background of the land-sea coexistence in the Bohai coastal region. Thirdly, compared to the other PBL parameterization schemes, the ACM2 scheme has the simulation results which are closest to real-time for the height of PBL, the profiles of potential temperature and water vapor mixing ratio. This is why the ACM2 scheme is more accurate in forecasting the heavy rainfall and the rain magnitude above it. Finally, the strength of the 700 hPa vertical movement simulated by the eight PBL schemes has an important effect on the northward typhoon precipitation, which basically determines the change trend of the hourly maximum precipitation and the relative magnitude of the regional average 24-hour cumulative precipitation.